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Rayleigh-Benard convection in horizontal layers of binary fluid mixtures heated from 
below with realistic horizontal boundary conditions is studied theoretically using multi- 
mode Galcrkin expansions. For positive separation ratios the main difference between 
the mixtures and pure fluids lies in the existence of stable three dimensional patterns 
near onset in a wide range of the parameter space. We evaluated the stationary solutions 
of roll, crossroll, and square convection and we determined the location of the stability 
boundaries for many parameter combinations thereby obtaining the Busse balloon for 
roll and square patterns. 



1. Introduction 

Convection in binary miscible fluids like ethanol- water, 3 He — 4 He, or various gas mix- 
tures shows a rich spectrum of pattern formation behaviour — see, e.g., (Platten & Legros 
1984; Cross & Hohenberg 1993; Liicke et al. 1998) for a review. The spatiotemporal prop- 
erties of convection in mixtures are more complex than those of one-component fluids due 
to the influence of Soret sustained concentration gradients. The structural dynamics of 
the concentration distribution in mixtures results from an interplay between three com- 
peting mechanisms: nonlinear advection and mixing, weak solutal diffusion, and the Soret 
effect. The latter generates and sustains concentration gradients in (linear) response to 
local temperature gradients. Without Soret coupling, i.e., for vanishing separation ratio 
%jj = any concentration fluctuation diffuses away. For ip ^ , however, the externally 
imposed vertical temperature difference across the fluid layer sustains via the Soret effect 
concentration variations against the action of advective mixing and diffusive dissipation. 

The concentration field changes the advective properties of mixtures via solutal buoy- 
ancy forces that enter into the momentum balance of the fluid. Thus, a concentration 
fluctuation directly influences the flow which in turn changes and mixes the concen- 
tration. In binary liquids, this nonlinear feed back is only weakly damped by diffusive 
homogenisation so that the concentration distribution shows anharmonic and boundary 
layer structures. Furthermore, it is ultimately this feed back that causes already right 
at onset convection patterns that cannot be seen there in pure fluids. Examples that 
occur depending on parameters are traveling waves of roll structures, standing wave os- 
cillations, and stationary squares. In addition mixtures show very interesting secondary 
structures close to onset: spatially localized traveling wave states, stationary crossrolls 
(CRs), and oscillations between squares and rolls or CRs can be seen for different pa- 
rameter combinations. Note that all of the aforementioned patterns are Soret induced by 
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the concentration field — they disappear in the pure fluid limit, tp — > 0, when switching 
off the Soret coupling to the concentration field. 

In this paper we are concerned with the case of a positive Soret effect, ip > 0, which 
causes the heavier (lighter) component of the mixture to be driven towards lower (higher) 
temperature regions. Therefore, heating a mixture with ip > from below establishes a 
stronger density gradient as in a pure fluid. The solutal contribution to the buoyancy 
increases the thermal destabilization of the fluid layer and convection starts at smaller 
temperature differences compared to a pure fluid. One commonly denotes the thermal 
driving regime with Rayleigh numbers R below the threshold R° c of pure fluids as the 
Soret regime and the regime above R° c as the Rayleigh region (Moses & Steinberg 1991). 
As a crude rule of thumb one can say that in the Soret region, R < square patterns 
are often observed in mixtures, whereas in the Rayleigh region stable rolls are found. 

There have been only few theoretical investigations aimed at explaining the transition 
scenario between squares at smaller R and rolls at larger R (Clune & Knobloch 1992; 
Muller & Liicke 1988). Recently we have elucidated this transition for a fixed wavenumber 
(Jung, Huke, & Liicke 1998). In this paper we compare the properties of square, roll, and 
CR patterns and we present a comprehensive linear stability analysis of rolls and a more 
restricted one for squares. We elucidate how the stability boundaries of rolls that have 
been determined by Busse and coworkers (Busse 1967; Bolton, Clever, & Busse 1985; 
Clever & Busse 1990) for pure fluid convection are modified by taking into account 
the influence of the concentration field in mixtures. We present for the first time a full 
numerical investigation of the stability behaviour of rolls and squares and present the 
stability balloons of these patterns in the k — r-plane for a wide range of fluid parameters. 

The paper is organized as follows: In § || we describe briefly the basics of convection in 
binary fluids and we explain the application of the Galerkin method to this particular 
system. In §p we describe and compare the stationary solutions for rolls, CRs, and 
squares. In §11 the Galerkin method is used for the linear stability analysis of roll and 
square patterns. We conclude in §|B| with a summary of our results. 



2. Mathematical foundations 

In this paper we investigate convection in horizontal binary fluid layers confined be- 
tween perfectly heat conducting, rigid, impermeable plates. Since the system and its 
basic equations are well known (Landau & Lifshitz 1966; Platten & Legros 1984), we 



summarize in § 2.1 only the necessary formulas for our investigation. Then we present 
relevant details related to the application of the Galerkin expansion technique to this 
system. 

2.1. System and basic equations 

We consider a horizontal layer of a binary fluid mixture of thickness d in a homogeneous 
gravitational field, g = —ge z . A vertical temperature gradient is imposed by fixing the 
temperature 

AT d . . 

T = T ±— at z = T - , (2.1) 

e.g., via highly conducting plates in experiments. Here we consider the plates to be 
infinitely extended, rigid, and impermeable. 

Convection is described in terms of the fields of velocity u = (u, v, u>), temperature T, 
mass concentration C of the lighter component, total mass density p, and pressure. In 
the balance equations connecting these fields we scale lengths and positions by d, time 
by the vertical thermal diffusion time d 2 / k, temperature by vn/agd 3 , concentration by 
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i/n/Pgd 3 , and pressure by pon 2 jd 2 . Here po is the mean density, k the thermal diffusivity, 
v the kinematic viscosity, and a and (3 are thermal and solutal expansion coefficients, 
respectively. Using the Oberbeck-Boussinesq approximation the balance equations read 
(Platten & Legros 1984; Hort, Linz, & Liicke 1992) 

V-u = (2.2a) 

(dt + u- V)u = -Vp + a [(8 + c) e 2 +V 2 u] (2.26) 

(dt+u- V)9 = Rw + V 2 9 (2.2c) 

(dt+u- V)c = Ripw + L (V 2 c-i/jV 2 9) . (2. 2d) 

Here 9, c, and p are the reduced deviations of temperature, concentration, and pressure, 
respectively, from the conductive profiles. 

The Lewis number L is the ratio of the concentration diffusivity D to the thermal 
diffusivity k, therefore measuring the velocity of concentration diffusion. The Prandtl 
number a is the ratio of the momentum diffusivity v and k: 

L=°-o-= V -. (2.3) 

The Rayleigh number R measures the thermal driving and the separation ratio ip mea- 
sures the strength of the Soret coupling between temperature and concentration fields 

R = <****, * = (2.4) 
vk a To 

Here To is the mean temperature and kr is the thermal diffusion ratio (Landau & Lif- 
shitz 1966). The driving forces entering into the momentum balance equation (2.26) 
are pressure gradients and the buoyancy caused by the temperature and concentration 
dependence of the density. 

The off-diagonal term -I i/; V 2 and the term Ripw in the concentration balance 
equation ( p.2c| ) describe the action of the Soret effect, i.e., the generation of concentration 
currents and concentration gradients by temperature variations. A Soret coupling tp > 
implies a positive Soret effect. In this case the lighter component of the mixture is driven 
into the direction of higher temperature thus increasing the density variations. 

The Dufour effect, i.e., the driving of temperature currents by concentration variations 
is of interest only in gas mixtures (Hort et aL1992). But even there it is often small (Liu 
& Ahlers 1997). 

2.2. Galerkin method 

To describe three-dimensional patterns with wavenumbers k x and k y each field X is 
expanded as 

X(x, y, z;t) = J2 X lmn (ty lk * x e imk yyf n (z) . (2.5) 

Ivan 

Here I and m are integers and the /„ form a complete system of functions that fits the 
specific boundary condition for the field X at the plates. To find suitable sets of functions 
/„ we introduce some new fields. First, two scalar fields $ and $ are defined via 

u = V x V x $e z + V x *e 2 . (2.6) 

The structures we want to discuss do not show a horizontal mean flow for mirror symme- 
try reasons. Then, ( |2.6| ) is the most general expression that fulfills the incompressibility 
condition ( 2.2cj ) (Clever & Busse 1989). The analysis of mean flow effects in the pertur- 



bations in discussed in section 
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Second, instead of c we use the field 

C = c-ip6 (2.7) 

that allows in a more convenient way to guarantee the impermeability of the horizontal 
boundaries: The diffusive part of the concentration current, driven by concentration gra- 
dients as well as by temperature gradients is given by — LV(c — rp0). At the impermeable 
plates the vertical component of this current vanishes which requires 

= d z (c - ipO) = d z Q at ^ = ±1/2. (2.8) 

The advective concentration current vanishes at the plates beca use there u = 0. The 
balance equation for £ is obtained by combining ( |2.2c| ) and ( 2.2c ). 
The boundary conditions for the fields 0, and ( read 

$ = <9 Z $ = * = 9 = d z ( = at z = ±1/2. (2.9) 

To expand the fields 6, £, and <& vertically we used different orthonormal sets f n (z) 
as follows 

T i/i <■ / \ f v / 2cos(n7rz) n odd , 

* and : f n (z) = i \- . } ' 2.10a 
|_ v2sm(n7rz) n even 

r 1 n = 

(■■fn(z)=\ V2 sin(mrz) n odd (2.106) 

[ a/2 cos(n7rz) n/0 even 

. „ f Cn+i (z) n odd 

* : = oT ■ (2-lOc) 

( n even 

Here C n and S n are Chandrasekhar functions (Chandrasekhar 1981). 
The balance equations for the new fields are 

<9 t A 2 # = ctV 2 A 2 * + {Vx [(u- V)u]} 2 (2.11a) 
<9 f V 2 A 2 $ = a {V 4 A 2 $ - A 2 [(1 + V) 6 + C]} - {Vx Vx [(u • V) u]} z (2.116) 
(d t + u- V)6» = -i?A 2 $ + V 2 (2.11c) 
(ft + u- V)C = iV 2 C-^V 2 . (2.11(f) 

Here A 2 = <9 2 +<9 2 . 

By inserting the ansatz ( |2.E| ) for each field into the balance equations and projecting 
them onto the basic functions one gets a nonlinear algebraic system of equations of the 
form 

A Kfi dfX^ — B Kfl Xfj, + C^vX^Xv . (2-12) 

For simplicity amplitudes are labelled here by a single Greek index and the summation 
convention is implied in ( [2. 12 ) with A KfJ _,B Kfl , and C KM „ being constant coefficients. 



The number of modes has to be truncated to get a finite number of equations as dis- 



cussed later on. For stationary convection structures the left hand side of (2.12) vanishes 
and the solution can be found using a multidimensional Newton method. 

2.3. Symmetries 

Symmetries of convective structures impose conditions on the fields and/or imply rela- 
tions between different modes of the fields thereby restricting the number of indepen- 
dent modes that are necessary to describe the patterns. For example, to describe two- 
dimensional roll patterns with k x = k, k y = that do not depend on y, all amplitudes 
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with 77i 7^ are set to zero in (2.5). On the other hand, square patterns are characterized 
by k x = k y = k and Xi mn — ±X m i n . But we also investigate three-dimensional CR 
patterns with k x — k y — k for which, however, Xi mn ^ ±X m i„. 

2.3.1. Stationary rolls 

To describe these two-dimensional structures one does not need the VP-field. Further- 
more, rolls are even in x with an appropriate choice of the plane x = 0. As a consequence 
of this mirror symmetry one has Xio n — AT_/on so that the lateral functions e ±lklx can 
be replaced by cos(klx). In addition the roll pattern is antisymmetric under reflection at 
the plane z = combined with a translation by half a wavelength in x-direction. This 
mirror glide symmetry enforces half of the amplitudes to be zero, e. g. all amplitudes 
$/on where / + n is an odd number. 

It is no accident that stationary roll patterns have these symmetries. That they are 
fulfilled at onset can be shown via a linear stability analysis of the conductive state 
(Hollinger & Liicke 1995). The subset of modes that obey these symmetries is closed 
in the sense that these modes do not drive others via nonlinear coupling. Thus, the 
observed roll solution remains symmetric as long as no symmetry breaking bifurcation 
occurs on the stationary roll branch with symmetry breaking modes becoming linearly 
unstable. Such instabilities are covered by our stability analysis. Moore, Weiss, & Willkins 
(1991) have discussed these symmetry breaking perturbations for free-slip and permeable 
boundary conditions. 

2.3.2. Stationary squares and crossrolls with k x = k y 

They have the same symmetry plane at x = as rolls and an additional mirror plane 
at y = 0. Furthermore, the squares and CRs show also a mirror glide symmetry. Here, 
however, the symmetry transformation consists of a reflection at the plane z = com- 
bined with a translation by half a wavelength in x- and y-direction. To describe these 
three-dimensional patterns the W-field cannot be neglected. We also mention that in 
contrast to the other fields 'J is odd in x and y and has positive parity under the mirror 
glide operation thereby reflecting the symmetries of the velocity field. 

For square patterns that are invariant under rotation by 90° in which the x- and 
y-directions are indistinguishable a further reduction of the number of mode occurs: 
amplitudes like $| ra „ and are the same. This is also true for 9 and £. Again ^ is 

different. Here * /m „ = 

2.4. Stability analysis 

To make a full stability analysis one has to check the stability of the patterns against 
perturbations with arbitrary wavevector de x + be y . To do so one has to introduce a 
Floquet term writing out the perturbation as 

SX(x, y, z; 1) = e idx + ib y e st £ 5X lmn e Uk " +imk y »/ n (z). (2.13) 

Irnn 

Such a perturbation is added to the known solution the stability of which is to be tested 
and inserted into the balance equations. After linearising and projecting one gets a linear 
eigenvalue problem of the form 

sA K „5X tl = B Kll 5X ll . (2.14) 

with constant coefficients A Kf + and B K/1 . The aforementioned solution, i.e., the convective 
structure described by it is stable if every eigenvalue s has a negative real part for every 
d and b. 
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The symmetry of the convective pattern discussed above can under some circumstances 
be used to get separated classes of possible eigenvectors representing the perturbations. 
That means the eigenvalue problem can be reduced to finding the eigenvalues of two 
matrices of about half of the size. Because evaluating the eigenvalues of a matrix is a 
0(iV 3 )-process this always implies a reduction of the computation time. 

2.4.1. Stationary rolls 

To perform the stability analysis of rolls one determines the growth behaviour of per- 
turbations of the form 

6X(x,y,z;t) = e ldx+ib » e st £ 5X l0 J lk * f n {z). (2.15) 

In 

Contrary to the roll solutions, their perturbations may contain a mean flow component 
that is discarded in (|2.6[ ). However, our perturbation ansatz ( [2.15| ) contains modes like 
5&ooie ldx+lbv Ci(z) describing a mean flow in the limit b,d 0, and the equation of 
motion for (J^ooi reduces to two independent equations for the mean flow as they are 
used e. g. in (Clever & Busse 1991). If there is a mean flow in the perturbations, then 
modes like the J^ooi-mode diverge when b, d goes to zero such that the long wavelength 
component of the velocity remains finite. However, in the parameter range explored here, 
the perturbations limiting the stability balloon of rolls have finite b or d and thus no mean 
flow. 

Because of the periodicity of the patterns in x and its mirror symmetry it suffices to 
consider d £ [0,k/2]. In y-direction, however, all perturbation wavenumbers, say, 6^0 
have to be investigated. 

The linear system of equations ( |2.15 ) always separates into two subsystems of per- 



turbations SX^ that belong to modes with amplitudes X^ that are antisymmetric (G- 
perturbations) or symmetric (G-perturbations) under the mirror glide operation (x, y, z) — 
[x + A/2, y + A/2, —z). E. g. all perturbations with amplitudes <5$;on with even I + n are 
G-perturbations, and the perturbations with odd I + n are G-perturbations. 

In general the perturbations ( [2.15| ) do not have a well defined symmetry under the 
mirror glide operation. This is only the case if d = nk. If n is even (odd) then the 
perturbations have the same (opposite) parity as the modes they belong to. Every G- 
perturbation can be written out as a G-perturbation (and vice versa) by choosing a new 
d 1 = d — k. Thus the distinction between G- and G-perturbations is well defined only 
for fixed d. 

It is possible to restrict the stability analysis to one set of perturbations by extending 
the investigated d-interval to [0, k]. Consider for example a G-perturbation with d £ 
[0, k/2]. It can as well be written as a G-perturbation with a new d' = d—k or equivalently 
d! = k — d. Thus one finds all G-perturbations with d 6 [0, k/2] again as G-perturbations 
with d! € [k/2,k]. We will therefore restrict ourselves to the set of perturbations that 
has negative parity under the mirror glide operation at d = 0, i. e., we investigate the 
G-perturbations in the whole interval [0, k]. 

In special cases the system of equations can be separated even further. For d = the 
perturbations can be divided into those that are symmetric and those that are antisym- 
metric under the operation x — > —x. Furthermore, if b = then the perturbations contain 
either no or only J^-amplitudes. 

2.4.2. Stationary squares and CRs with k x = k y 

Because the amount of computational power needed to make a full stability analysis 
of these three-dimensional structures is too large, we will discuss perturbations only for 
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periodic boundary conditions, i. e., d = b = 0. Here again a separation of perturbations 
is possible into those that change sign or not under the mirror glide operation. 

Furthermore, the stability problem is invariant under x — > —x and y — > —y. Thus one 
can distinguish between perturbations that are even in x and y, odd in x and y, or even 
in x and odd in y (or equivalently odd in x and even in y). If the perturbations have 
the same symmetry in both directions one can in the case of squares finally make use 
of a last symmetry property and separate between perturbations that are symmetric or 
antisymmetric under the exchange of the x- and y-dircction. This is not possible for CR 
patterns because of their lack of the x <-> y symmetry. 

As we will see in section the destabilising perturbations of squares fall into the 
subclass that is even in x and y, and therefore do not drive a mean flow. 



3. Properties of the patterns 

Squares, CRs, and rolls are realized as stable convection structures somewhere in pa- 
rameter space. Furthermore, if L is sufficiently small these three patterns appear for 
fixed L, a, if), and A; as global attractors at different Rayleigh numbers (however, CRs 
coexist bistably with oscillations in a small i?-interval (Jung et aZ.1998)). So it is easy 
to enforce their stable experimental realization successively by increasing the Rayleigh 
number beyond the onset of convection. 

To understand the behaviour of the patterns it is useful to consider the driving region 
near onset and the Rayleigh region separately. In the pure fluid the critical point is at 
a Rayleigh number R° = R c {ip = 0) = 1707.762 and a wavenumber fc° = k c (i]j = 0) = 
3.117 (Chandrasekhar 1981). In binary mixtures with positive separation ratios that we 
are dealing with, the critical Rayleigh number is smaller, R c (ip > 0) < R®, since the 
solutal contribution to the quiescent state's buoyancy force enhances the latter. Thus 
a smaller thermal driving, i.e., a smaller Rayleigh number suffices to reach the critical 
buoyancy force size for onset of convection. The critical wavenumber is also somewhat 
lower: k c (ip > 0) < fc° (Knobloch k Moore 1988). 

When presenting our results we shall use the reduced Rayleigh number 

r = l < 31 > 



and the reduced distance 



from threshold r c = jjj; 



R r , s 



3.1. Small- amplitude convection — amplitude equation 

Very close to onset the wavenumber dependent bifurcation behaviour of rolls and squares 
can in pure fluids as well as in binary mixtures be described most simply by two coupled 
cubic amplitude equations of the form 

To d t A = eA + e o (d x -^d 2 )j A-\A\ 2 A-f\B\ 2 A (3.3a) 
r d t B = eB + e (dy - B - \B\ 2 B - f\A\ 2 B . (3.36) 
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Clune and Knobloch used such equations without the spatial derivative term (Clune & 
Knobloch 1991). For a review of such amplitude equations and how they are related to 
the basic equations see, e. g., (Cross & Hohenberg 1993). 

The amplitudes A and B of the eigenfunctions of the linearized hydrodynamic field 
equations correspond quite well to the leading amplitudes 

wioi = & 2 $ioi = k 2 J §(x, y, z) e ikx C\ (z) dxdydz (3.4a) 
w 011 = /c 2 $on = k 2 J §(x, y, z) e iky d [z) dxdydz . (3.46) 

of the vertical velocity field or the corresponding amplitudes of the temperature field. 

The two roll solutions of (3.3) are A = \A^\e 1 ^ k ~ kc " >x , B = for rolls with wavevector 
k = ke x and A = 0,B = |_B R |e'( fc ~ fcc * )a for rolls with wavevector k = ke y . Squares are 
described by the solution A = \ A s \e^ k ~ k ^ x , B = \B s \e i( - k - k ^ y with A s = B s . Crossrolls, 
i. e. solutions with finite \A\ ^ \B\ do not exist in (3.3). 

If / > — 1, the square solution exists besides the roll solution for e above the neutral 
stability curve 

estab(k) =£(k~k c ) 2 . (3.5) 

We will always consider roll patterns with Br = and Ajj ^ 0. For the roll solutions 
one has 

\A R \ 2 = e-e stab (k) , (3.6) 

and for squares one finds 

\A S \ 2 = \B S \ 2 = IZiggp ( 3. 7) 

so that 

|^| 2 = (1 + /)|A S | 2 . (3.8) 

A stability analysis of the roll and square solutions shows that squares (rolls) are stable 
(unstable) whenever |A^| 2 < 2|As| 2 , i.e., for — 1 < / < 1 and they are unstable (stable) 
whenever |Ar| 2 > 2|v4 s | 2 . 

3.2. Full Galerkin expansion 

Both square and roll solutions exist for every r above onset r c (ip,L) and the mode 
intensity |u>ioi| 2 for rolls is always greater than for squares. Hence the parameter / in 
the amplitude equation has to be greater than according to (3.8). In figure [l] we present 



a plot of the leading contribution |wioi| 2 + |won| 2 to the vertical flow intensity w 2 versus 
r as obtained from the full Galerkin expansion for several parameters. 

By comparing our results from the full Galerkin expansion with the amplitude equation 
approximation we verified that the latter works well near onset r c (f/ ; ) which in figure 
|l| lies significantly below 1 - outside the plot range of figure [l|. However, at larger r 
when approaching the Rayleigh regime the amplitude approximation loses its validity: 
\A\ 2 + \B\ 2 continues to grow linearly with r with its initial slope at onset whereas 
|it?ioi | 2 + l^on | 2 strongly curves upwards in the Rayleigh regime. Furthermore, while 
l w ioi| 2 < l^foil 2 + l w onl 2 = 2|wf i| 2 holds close to onset — as required within the 
amplitude equation approximation for squares to be stable — the full solutions at larger 
r are such that Iwf^l 2 > |wf i| 2 + l^onl 2 without, however, rolls becoming stable, what 
happens at even higher r. 

In the Rayleigh region the amplitudes become comparable with those of the pure 
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fluid. This is because the concentration field gets more and more advectively mixed and 
equilibrated and loses therefore its influence on the convection. Pure fluid convection, 
ip = 0, can be described close to onset r c (ip — 0) = 1 by a cubic amplitude equation. 
However, the initial slope of w 2 versus r is much greater than for the binary mixtures 
shown in figure ^ with ip > 0. The transition between Soret and Rayleigh region is 
especially sharp at small L - c.f. the bifurcation diagram for L = 0.01 in figure 0(b). We 
will observe such a behaviour again when we consider the stability of the patterns. 

In the r-range investigated in this paper the field amplitudes ^loi and $on that 
are kept in the amplitude equation approximation are also the leading ones in the full 
Galerkin expansion. To describe only the fixed point solutions it would indeed be sufficient 
to approximate the velocity field of rolls (squares and CRs) by one mode (two modes) 
only. But this restriction does not suffice when the linear stability is investigated. 

Although more than two amplitudes are needed, a good representation is easier to 
achieve for the velocity field than for the £ _ and 0-fields. For them very many modes 
are needed if L is small and r is large. This is a region of the parameter space where 
the concentration field shows narrow boundary layer behaviour which has to be resolved 
properly. In addition a consistent description of the temperature field requires then also 
— despite the fact that it is rather smooth — high S-modes as discussed by (Hollinger 
& Liicke 1998; Hollinger 1996). 

We followed (Clever & Busse 1989) when defining our truncation prescription for the 
Galerkin expansion. We defined a maximal mode index N and neglected all modes Xi mn 
with \l\ + \m\ + n > N. We took the smooth behaviour of the velocity field into account 
by defining two different indices: N\ for the $- and ^-fields and N2 = 2Ni for the 0- 
and (-fields. For the most anharmonic roll structures at r « 1.5, L < 0.01, and ip = 0.15 
that we have investigated expansions up to N2 = 40 were needed. This is much more 
than for pure fluids, where truncations with N 8 are sufficient to describe the stability 
behaviour quantitatively even at large r. Since the structure of squares is somewhat 



smoother than the one of rolls (cf. § 3.4) and since they exist stably only at small r in 
or near the Soret region, we needed only N2 ^ 20 for squares. Also the CR structures 
could be described well with such a truncation near the first bifurcation point at small 
r, where they behave square-like. But in order to resolve the CR structures also close 
to the second bifurcation point at larger r where the CR solution merges into the roll 



solution (cf. § |3.3| ) more modes would have been necessary. 



3.3. Bifurcation behaviour of rolls, squares, and crossrolls 

Figure^ shows a typical bifurcation diagram for a parameter combination where the three 
stationary patterns can be found. It also contains information on the stability of these 
patterns. Crossrolls exist only in a finite r-interval. The CR solution branches emerge 
out of the square branch slightly above r — 1. At r ^ 1.36 the CR solution disappears 
when, e.g., the amplitude (downwards pointing triangles in figure ||) becomes zero 
and the CR branch for (upwards pointing triangles) ends on the roll solution branch 

w ?oi- 

On the other hand, roll as well as square solutions exist for all r ^ r c of figure g. 
Initially at onset the latter are stable and the former are unstable. For the parameters of 
figure H squares lose their stability in a Hopf bifurcation at r ~ 1.11 to oscillations which 
eventually undergo with increasing r a subharmonic bifurcation cascade that is termi- 
nated when the CR states have become sufficiently attractive to quench the oscillations. 
For other parameter combinations, in particular for larger L there are no oscillations and 
the squares transfer their stability directly to CRs (Jung et al. 1998). 
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3.4. Structural properties of roll and square fields 

In figure ^ we show the concentration distribution of square convection for two parame- 
ter combinations that are representative for liquid and gas mixtures. This plot and the 
concentration field structure of rolls and squares in a vertical cross section shows a char- 
acteristic boundary layer and plume behaviour at small L. Such structures occur when 
advective mixing is large compared to diffusion in the bulk of the fluid. Consequently the 
boundary layers and plumes are more pronounced in rolls than in squares since w 2 R > w| 



as discussed in § 3.2. Thus squares with their broader boundary layers are much smoother 
structures than rolls at the same parameters. 

The practically harmonic velocity and temperature fields are not shown. For squares 
they resemble the fields of a linear superposition of two perpendicular sets of rolls. 

The Nusselt number TV is roughly the same for rolls and squares. Close to onset 

TV -lex Koi| 2 + Kn| 2 , (3-9) 
and the stable structure h as th e higher Nusselt number there, i.e., Nr < N$ in accordance 



with the inequalities of § |3^1| predicted by the amplitude equation. Further away from 



onset, however, one has N$ < Nr thus reflecting the magnitude relations of w 2 discussed 



in §y. 

In figure || we show the mixing parameter 



M f*P ,3.10) 



co i id 



It is defined by the mean square of the deviation, SC — C — Co, of the concentration 
from the spatial mean, Co = (C), reduced by the concentration variance in the quiescent 
conductive state. Note that M is nearly the same for the very different concentration 
fields of rolls and squares, if L is not too large (figure 0). 



4. Linear stability analysis of rolls and squares 

4.1. Instability mechanisms of rolls 

The stability boundaries of roll patterns in pure fluids are known since the pioneering 
work of (Busse 1978). At small Rayleigh numbers there exist five different instability 
mechanisms giving rise to five different stability boundaries that limit the region of sta- 
ble rolls in the (R, k, c)-parameter space. At small Prandtl numbers the Eckhaus (EC), 
the skewed varicose (SV), and the oscillatory mechanism (OS) are the important insta- 
bilities (Clever & Busse 1990). At higher Prandtl numbers the zigzag (ZZ) and the CR 
mechanisms dominate (Busse 1967). Properties and symmetries of these perturbations 
are discussed in (Bolton et al. 1985). All these five instabilities of roll patterns can also 
be found in binary mixtures. 

In pure fluids and binary mixtures there exist always perturbations of the form (2.15) 
with b = d = against which a roll pattern is only marginally stable. Such a perturbation 
has no ^/-component and is odd in x. It reflects just an infinitesimal shift of the whole 
pattern in the ^-direction. Therefore such a perturbation has an eigenvalue s — 0. For 
later discussion we point out here that this particular eigenvalue is connected to nearly 
all instabilities: the perturbations causing them and the associated eigenvalues evolve 
smoothly into the lateral shift when one moves in the d— b plane from the d— b coordinates 
that locate the instability to the origin d = b = 0. 



In the remainder of § 4.1 we will briefly characterize the properties of the aforemen- 



tioned five perturbations before presenting our results of the stability analyses for rolls 
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in §4.2 and for squares in §4.3. We begin with the three types (EC, ZZ, and CR) 



that touch the critical point e = 0, k = k c and that exist also within the two coupled 



amplitude equations (3.3a - 3.3 b) 



4.1.1. Eckhaus instability 

Perturbations of the Eckhaus type are most critical at b = and have no ^-component. 
Thus they can be described as purely two-dimensional. The EC instability tends to es- 
tablish a new set of rolls with a better wavenumber in the direction of the wave vector 
of the unstable roll pattern. Within the amplitude equations EC perturbations can be 
represented as variations of the A-amplitude that depend only on x. Instability occurs 
here for e and k such that e sta b(k) < e < eEc(k) = 3e sta b{k). At d — ► the EC perturba- 
tions reduce both in the amplitude equation approximation and in the full hydrodynamic 
equations to the lateral shift. For symmetry reasons s ~ d 2 near d = 0. In the case of 
EC instability (stability) there is a minimum (maximum) of s at d = 0. It is therefore 
sufficient to investigate the pattern at a single point on the d-axis near d = numerically 
in order to determine the stability behaviour against EC perturbations. But to find the 
most critical value of d an evaluation and interpolation of s along the d-axis is necessary. 

4.1.2. Zigzag instability 

Zigzag perturbations have d = and fall into the subclass of perturbations that are 
odd in x. In the amplitude equations they show up as y-dependent perturbations in 
the ^4-amplitude. They cause the growth of a new set of rolls that has always a greater 
wavenumber than the original set. Consequently they confine the region of stable rolls 
on the small-fc side. The amplitude equations predict a ZZ instability for all k < k c . 
Like EC perturbations the ZZ instability reduces to the lateral shift when 6-^0 and the 
question of stability can be answered at a single point near d = b = 0. 

4.1.3. Crossroll instability 

It occurs when roll-like perturbations perpendicular to the existing pattern can grow. 
In the amplitude equations ( |3.3q - 3.36) they are described as perturbations in the 



amplitude B when A describes the stationary roll pattern. Rolls are CR-unstable for 
£stab(k) < e < e C R(k) =//(/- 1) e sta b(k) when / > 1. In this case e C R(k) can be above 
or below the Eckhaus boundary CEc{k) — 3e sta b(k) depending on whether / < 3/2 or 



not. However, if / < 1 (which is the case when squares arc stable c. f. §3.1) then 
rolls are Ci?-unstable for all e > e s tab{k) within the amplitude equations approximation 
( [3.3a - 3.36). Note that the case / < 1 does not occur in pure fluids but it can occur in 



binary mixtures when squares are stable at onset. 

In the multi-mode Galerkin expansion the leading mode in the velocity field of the CR 
perturbation has the form 

8^ 00l e ihv C 1 {z) , (4.1) 
where near the critical point b ss k c . This is a mode of the G class of perturbations that 



are symmetric under the mirror glide operation as discussed in § 2.4.1. Being a member 
of the G class this perturbation cannot be connected smoothly in the d — b plane to 
EC or ZZ perturbations since the latter belong to the G class of perturbations that 
are antisymmetric under the mirror glide operation. Since we decided to transform all 



perturbations into the G class by a shift d' = d — k as explained in § 2.4.1 we have to 



rewrite the above CR perturbation (4.1) into the form 



SQ-wie^-^+^diz) . (4.2) 



Writing out the CR perturbations in this form one finds again that the corresponding 
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eigenvalue is connected to the eigenvalue of the lateral shift, only taken at a different 
d = k. 

Since the CR instability does not occur at arbitrarily small b one therefore has to test 
several values of b and then apply an interpolation procedure. 

4.1.4. Skewed varicose instability 

This instability is not captured by the simple amplitude equations. The SV boundary 
confines the stability balloon on the large-/c side. When crossing this boundary the per- 
turbation tends to replace the original set of rolls by a new set with smaller wavenumbcr. 
The eigenvalue has its maximum at d ^ and 6^0. The SV instability too reduces to 
a lateral shift at d, b — > and can be found at infinitesimal small b and d. To find the 
SV stability boundary of rolls one has to find the maximum of the eigenvalue on a line 
between the d— and 6-axis near the origin. 

4.1.5. Oscillatory instability 

The oscillatory instability is the only instability with complex eigenvalues at small 
r. The perturbation is most critical at d — and finite b. It is even in x. The pair of 
complex conjugate eigenvalues undergoes near b — a collision and generates two real 
eigenvalues. One of these stationary perturbations transforms into the lateral shift. The 
real eigenvalues are negative (if the pattern is stable against ZZ perturbations), so the 
search for the OS instability requires an evaluation and interpolation along the fr-axis as 
for the CR instability. 

4.2. Stability boundaries of rolls 

We saw that the EC, ZZ, CR, SV, and one eigenvalue of the OS instability are connected 
to the lateral shift at b = d = 0. Figure || shows an example of the most dangerous 
eigenvalue in the (d, &)-plane. Along the d-axis the value of s goes down from s = at 
the origin. The pattern is therefore EC stable. The maxima at d = and d — k show 
that it is unstable against ZZ and CR perturbations. An oscillatory instability does not 
occur here for these parameters since the eigenvalue is always real. There is also no SV 
instability. The latter would cause a relative maximum between the d— and &-axis. 

Besides the perturbation at d = b = which describes a lateral shift in all fields there 
is another important location in wavenumber space that is of relevance for concentration 
field perturbations and which thus is specific to binary mixtures. It lies at d — k, b = 
and describes a change of the mean concentration. This is most easily understood by 
transforming this perturbation from the G class into the G class where it then occurs at 
d = b = 0. Here it consists of one single mode <5(ooo- This mode is constant in the fluid 
layer and describes a change in the average concentration. But since only derivatives of the 
(■-field show up in the balance equations such a mode has no influence so that s = 0. The 
reason why this mode has to be taken into account is that for d, b ^ it describes a long- 
wavelength perturbation 5CQQ0 etdx+lby that is indeed of physical importance. Because 
the average concentration is fixed, a divergence in the <5(ooo -am phtude (as in <5<I>ooi 
for mean flow) cannot occur when approaching b, d = 0. While the eigenvalue of the 
CR perturbation is connected to the above described zero eigenvalue the CR instability 
occurs always at finite b. Only far in the unstable region the CR eigenvalue becomes 
positive at arbitrary small b. However, we found no roll instability to occur directly near 
this point whereas EC, ZZ, and SV instabilities are realized near the origin d = b = as 
modifications of the lateral shift. 

Figure || shows the CR instability boundaries of rolls in the full Galerkin model in a 
parameter interval where an exchange of stability between rolls and squares at onset is 
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predicted in (Clune & Knobloch 1991). One sees that the curvature of the CR boundary 
at the critical point diverges — as is also predicted by the amplitude equation (3.3c - 3.3&) 
for / = 1 — when this exchange occurs. For the parameters a = 10, ip = 0.01 of figure [j] 
the exchange occurs at L = 0.2. Decreasing L below this value the neutral stability curve 
(dashed line in figure ^ drops further down in r (not shown in figure ^|b) while the CR 
instability boundary moves up in r. The r-range between these two curves locates stable 
squares. In such a situation where squares are stable at onset the amplitude equations 
predict rolls to be CR-unstable everywhere while in the full equations rolls become stable 
against CR-perturbations above the solid lines in figure |^. 

The rolls could still be unstable there against other perturbations but we found the 
minimum of the CR boundary always to be the minimal Rayleigh number for stable 
rolls to exist. To know the location, r c CR , of this minimum of the CR boundary for 
several parameters is therefore of interest. We have determined r c CR and the associated 
wavenumber k c CR of the most dangerous CR perturbation for ip = 0.01 and ip = 0.15 and 
presented the results in figures [?] and 0, respectively. They show that rolls are CR unstable 
at onset, r c (dashed lined in figures 0a and ||a), for small L and large a. However, for 
larger L being typical for gas mixtures, rolls are stable at the critical point. 

At iji = 0.15 the minimum r c CR of the CR boundary and its wavenumber location k c CR 
strongly increase for small a and L — cf. the small- L variation of the respective curves 
at a = 0.2 in figures |^a,b. r c CR seems to diverge here, but a detailed inspection shows 
that this is not true. The two branches of the CR boundary meet again at higher r 
near the apparent divergence. They limit in the k — r plane an oval region of CR-stable 
rolls from below and also above. (The upper part of this region can be unstable against 
OS-perturbations, though.) By reducing L r c CR remains finite but the oval region gets 
smaller until the region of stable rolls vanishes. The experimental observation of this 
behaviour might be difficult because it occurs in a region of the parameter space that is 
not accessible by ordinary fluid mixtures. 

We did also investigate the behaviour of b^ R , the wavenumber of the critical pertur- 
bation at (kQ R ,rQ R ). Only in the Rayleigh region the values of k^ R and b c CR are near 
k c = 3.117. But in general k c CR ^ b c CR . The linear analysis gives thus a hint for the 
existence of patterns with k x ^ k y in the Soret region. 

We have calculated all stability boundaries at small r for different values of L, ip, and 
a. Concerning the stability behaviour of the roll structures one sees that in the Rayleigh 
region, r > 1, where the concentration field is nearly uniform the binary mixture behaves 
like a pure fluid. And the transition between Soret and Rayleigh region is very sharp 
at small L and ip. An example for such a behaviour is given in figure ^[ Here only the 
EC, CR, and ZZ boundaries are of importance. In the Rayleigh region r > 1 of figure || 
the CR, ZZ, and EC boundaries of the mixture (full lines with circles, triangles, and 
squares, respectively) are lying close to the corresponding boundaries of the pure fluid 
(long-dashed lines). Note in particular the vase-like form of the EC boundary rEc(k) and 
the dent in the ZZ boundary rzz{k): Close to onset (r c ~ 0.6, k c ~ 2.6 in figure ^), i.e., 
in the Soret regime rEc(k) opens up parabolically and rzz(k) comes out of the critical 
point linearly with positive slope. However, in the crossover range r ~ 1 between Soret 
and Rayleigh regime the curve rEc{k) pinches inwards and develops a waist such as to 
follow in the Rayleigh regime the parabolic shape of the EC curve of the pure fluid that 
starts out at ~ 3.1, = 1. Similarly rzz{k) bends in the crossover range towards 
small k to follow then the ZZ boundary of the pure fluid that shows negative slope. 

Figures |K)|-|l2| show that this sharp transition between Soret and Rayleigh region that 
causes the vase-like structure of the EC boundary and the sharp bend of the ZZ boundary 
does not occur when L is greater. For L = 0.01 and ip > 0.08 the critical point lies at 
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k = 0. In this case the amplitude equations are not applicable. We found the ZZ boundary 
and the left EC branch to go to r — oo at small k. But the right EC branch does still 
meet the critical point. 

At low Prandtl numbers the roll solutions have a finite region of stability but their basin 
of the attraction seems to be very small since typically spiral defect chaos is observed here 
in experiments (Liu & Ahlers 1996). But it was shown that rolls could also be observed 
with special experimental procedures (Cakmur et al. 1997). 

We tried to find laws that connect the position of the boundaries and the fixed point 
solutions. We found no easy connection between these values for CR, ZZ and SV. But 
the OS boundary for fixed a seems to be independent of L and ip in the plane of k and 
the convection amplitude u>ioi instead of in the k—r plane (figure |l3| ). 

For the EC boundary a more complicated procedure is needed. We define an effective 
control parameter by linear interpolation of the values for the convection amplitudes 

^eff(r)=wl Q1 , (4.3) 

and an effective reduced Rayleigh number 

r*ff =£eff+r° stab (k) . (4.4) 

Within the amplitude equations it is just e e // = e. Plotting the EC boundary in the 
k — r e ff plane instead of in the k — r plane shows in the Rayleigh region only a dependence 
on a but not on ip and L, c.f. figure |lj. However, this procedure does not hold in the 
Soret region. 

4.3. Stability boundaries of squares 

Performing the stability analysis of squares we had to restrict ourselves to the case 
d = b = where the perturbations separate into different symmetry classes. An analyti- 
cal investigation of long wavelength perturbations of squares near the critical point can 
be found in (Hoyle 1993). Because both squares and rolls can be described as even in 
x and y and as mirror glide antisymmetric, one expects a perturbation that destabilises 
the squares and favours the rolls to fulfill these symmetries, too. However such a pertur- 
bation should break the symmetry x <-> y. We actually always found the most critical 
perturbation to fall into this symmetry class. Other perturbations that break the mirror 
symmetry in x— or y-direction are less critical. 

Figure |l5| shows typical examples for the stability region of squares. The left and 
right side of the stability boundaries should not be taken too serious — presumably 
square structures are destabilised earlier by instabilities with finite b or d that tune the 
wavenumber and that are not considered here. 

Even if squares are stable at onset they always lose their stability against a roll pattern 
at higher r. Furthermore, for certain parameters there does also exist a band of i — values 
where neither squares nor rolls are stable. Within this band three-dimensional CRs that 
break the x <-> y symmetry can be stable. For certain parameters we found stationary 
CRs only. For others oscillating CRs structures appeared as well. We have studied these 
structures in more detail in (Jung, Huke & Liicke 1998). 

5. Conclusion 

We investigated roll, crossroll, and square convection in binary mixtures for a wide 
range of parameter combinations using a multi-mode Galerkin method. All these pat- 
terns are realized as stable convection structures somewhere in parameter space. The 
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bifurcation behaviour of rolls and squares can be modelled near onset by amplitude 
equations, which, however lose their applicability in the Rayleigh region. Moreover, the 
CR solution that connects the roll and the square branch at Rayleigh numbers r w 1 is 
absent in a simple ansatz of cubic amplitude equations. 

We compared the Nusselt number and also the mixing parameter of the full numerical 
solutions for rolls and squares. We found these global properties of the convective states 
to be approximatively the same for these patterns despite the qualitatively different 
structures of the concentration fields of squares and rolls: Rolls show a strong boundary 
layer behaviour at small L and high r, so that compared to pure fluids much more 
modes are needed to describe them. Squares, on the other hand, show a less pronounced 
boundary layer behaviour making it easier to determine these three-dimensional solutions 
numerically 

In the main part of the paper we investigated the linear stability of rolls and squares. 
To that end we performed a full and unrestricted stability analysis for rolls using arbi- 
trary perturbations and in addition a stability analysis of squares that uses the periodic 
boundary conditions of the squares also for the perturbations. The result is that squares 
are stable in the Soret region if L is sufficiently small. But they always lose their stability 
at higher r. Rolls are stable at onset at higher L and in the Rayleigh region where squares 
are unstable. Typically however, there is a finite interval of r- values in between where 
neither rolls nor squares are stable. In this interval cither stationary or oscillatory CR 
patterns are observed. 

The analysis of the rolls shows that in the explored parameter range only the basic 
mechanisms of instability occur that are already known from the pure fluid, namely 
the Eckhaus, zigzag, crossroll, oscillatory, and skewed varicose mechanism. When the 
Soret region is small, the stability balloon of the mixtures resembles the Busse balloon 
for the pure fluid. However, at small L when the Soret region is large the situation is 
different. The fixed point solutions show a sharp transition between the two regimes. 
The convection amplitudes are very small in the Soret regime. But near r = 1 they 
increase strongly and e. g. the Nusselt number becomes comparable to that of the pure 
fluid convection. The stability boundaries of roll convection show a similar transition 
here. In the Rayleigh region the boundaries are close to the boundaries of the pure fluid. 
But upon reducing the Rayleigh number the boundaries begin to deviate from their pure 
fluid counterparts. Near onset they finally agree with the predictions of the amplitude 
equations for the mixtures. 

The EC boundary gives a typical example for this behaviour. For high L it has just a 
normal parabolic shape. But for small L one observes a qualitatively different behaviour. 
When approaching r = 1 from above the two branches do not meet near r = 1 as in 
the pure fluid but begin to separate again until r gets small enough for the amplitude 
equations to become valid. A typical vase-like shape results. 

A detailed inspection of the point where the rolls become unstable at onset shows that 
the CR boundary is responsible for this loss of stability. As long as the rolls are stable 
at onset, the CR boundary touches the critical point, as predicted by the amplitude 
equations. At a certain point the amplitude equations show a global loss of stability 
against CR perturbations. In the full Galerkin analysis this loss of stability does not 
occur. The CR boundary disconnects from the neutral curve and the rolls remain stable 
not at onset but at higher r. 

Both the amplitude equations and the full Galerkin expansion show that square struc- 
tures get stable at onset when roll structures lose stability there and vice versa. The 
stability domain of squares lies mainly in the Soret region and it is separated from the 
region of stable rolls at higher r by the region of CR structures. The fact that square 
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patterns lose their stability already at relatively small r was another favourable property 
which delimited the requirements for a numerical analysis of these patterns. 

The perturbations against which square solutions are unstable show those symmetries 
that squares, rolls, and CRs have in common. But they break the x <-» y symmetry, 
which exists only for squares. The resulting stability boundary delimits the region of 
stable squares in the k-r plane not only on top but also at the sides at smaller and larger 
wave numbers. It should be noted that our analysis of squares does not cover instabilities 
that retune the actual wavenumbcr of the considered pattern and which might occur 
earlier than the wave number preserving instabilities investigated here. 

This work was supported by the Deutsche Forschungsgcmeinschaft. 
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Figure 1. Flow intensity of stationary squares and rolls. Shown is the contribution 
fc 2 (|'I , i()i 2 + <l?oii | 2 ) from the leading modes to w 2 versus reduced Rayleigh number. Parameters 
are a = 1, L = 0.01 (a), and ip = 0.08 (b). The thin curves denote the roll solution in pure fluids 

y> = o). 



Figure 2. Representative bifurcation diagram of stationary patterns. The two most important 
modes of the velocity field are plotted versus r. Squares, circles, and triangles denote the square, 
roll, and CR solutions, respectively. Upwards (downwards) pointing triangles refer to, say, |wuu | 2 
(l w on| 2 ) f° r CRs with dominant 101 mode. Closed (open) symbols denote stable (unstable) 
structures. In the interval 1.11 < r < 1.15 the stationary solutions are unstable but an oscillatory 
solution is stable (Jung et aZ.1998). Parameters are L = 0.01, a = 10, ij) = 0.15, and k = k°. To 
get a consistent bifurcation diagram all solutions had to be calculated with the same truncation 
level N2 = 16 (see text). However, this is not sufficient to describe rolls quantitatively. 



Figure 3. Structural properties of square (a-d) and roll (e,f) convection for representative liquid 
parameters (L — 0.01, a = 10, left column) and gas parameters (L = a = 1, right column) at 
r — — 0.15. In (a,b) the concentration distribution of squares at mid height, z = 0, is shown. 
In (c-f) we show the concentration distribution in a vertical cross section at y = 0. The largest 
vertical upflow is at x = y — 0. 



Figure 4. Mixing parameter M (3.1C) versus r for rolls (circles) and squares (squares). 



Parameters are k = 3.117, a = 1, and L = 0.1(a), ij) = 0.08(b). 



Figure 5. Variation of a representative eigenvalue s as a function of perturbation wavenumbers 
d, b. The positive maximum at d — and b ~ 1.5 (d = k and b « 3.0) implies ZZ (CR) instability. 
Parameters are r = 1, k = 2.7, a — 10, ip — 0.08, and L — 0.1. 



Figure 6. CR instability boundaries (solid lines) of rolls in the Soret region for a = 10, ip = 0.01 
and several values of L. Rolls are stable against CR-perturbations above the solid lines, (a): For 
L ^ 0.2 the CR boundaries touch the neutral curve (dahed lines) in the critical point, (b): For 
L < 0.2 the neutral curve goes further down (not shown) and is disconnected from the CR 
boundaries. Then rolls are not stable at the critical point anymore. 



Figure 7. (a): The minimum of the CR-boundary r c CR denoting the smallest r where CR-stable 
rolls exist as a function of L for several values of a and V = 0.01. (b): The corresponding 
wavenumbers kQ R . 



Stability boundaries at positive separation ratio 
Figure 8. Same as figure m, but for ip — 0.15. 
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Figure 9. Stability boundaries of rolls. Solid lines with open symbols refer to a mixture with 
a = 7, ip — 0.01, and L — 0.025. The corresponding boundaries in a pure fluid with a — 7 are 
the long dashed curves. Dotted curves labelled EC(AE) and ZZ(AE) are the predictions of the 
cubic amplitude equations for the EC and ZZ boundaries of the mixture, respectively. 



Figure 10. Crossections of the stability balloon of rolls in the k — r plane at a = 0.1. Open 
circles: CR, open squares: EC, open triangles: ZZ, filled circles: SV, filled squares: OS. S denotes 
the region of stable rolls. For tp — 0.15, L = 0.01 there are no stable rolls. 



Figure 11. Crossections of the stability balloon of rolls in the k — r plane at a — 1. Symbols 

have the same meaning as in figure ^ 



Figure 12. Crossections of the stability balloon of rolls in the k — r plane at a = 10. Symbols 

have the same meaning as in figure ^ 



Figure 13. (a): OS boundaries ros{k) for the nine parameter combinations of figure [To| versus 
k. (b): Amplitude square |wioi| 2 (fc) of the leading velocity mode of the roll patterns at the OS 
boundaries shown in (a). Thick lines in (a) and (b) refer to the pure fluid at the same a — 0.1. 



Figure 14. (a): EC boundaries rEc(k) for the 27 parameter combinatio ns of fig ures O - [l^ 



versus k. (b): EC boundaries of (a) in the k — r e jf plane defined by (4.3, 4.4). Only data 
rEc{k) > r° atab {k) in the Rayleigh region are shown. The location of the EC boundary depends 
only slightly on a for a > 1. Therefore, only two different boundaries appear in (b). The inset 
shows the definition of e e //- 



Figure 15. Crossections of the stability balloon of squares in the k~r plane obtained at a — 0.1 
(open squares) and a = 10 (open circles) from a restricted stability analysis as explained in the 
text. S denotes the region of stable squares. For a = 0.1, L — 0.1, and ip = 0.01 there are no 
stable squares. 
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